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We study the quantum dynamics of Gaussian wave packets on star graphs whose arms feature 
each a periodic potential and an external time-dependent field. Assuming that the potentials and 
the field can be manipulated separately for each arm of the star, we show that it is possible to 
manipulate the direction of the motion of a Gaussian wave packet through the bifurcation point by 
a suitable choice of the parameters of the external fields. In doing so, one can achieve a transmission 
of the wave packet into the desired arm with nearly 70% while also keeping the shape of the wave 
packet approximately intact. Since a star graph is the simplest element of many other complex 
graphs, the obtained results can be considered as the first step to wave packet manipulations on 
complex networks. 


I. INTRODUCTION 

Recent experimental advances have led to very pre¬ 
cise manipulations of atoms in optical lattices [1] and of 
wave packets in waveguide arrays [2, 3]. Fundamental ef¬ 
fects such as Bloch oscillations for wave packets in tilted 
lattices have been experimentally confirmed by several 
groups [4-7]. Moreover, it has been shown that coher¬ 
ent control over the wave packets is possible by varying 
the external field: Such variations can be discrete [8, 9] 
or continuous [10-13]. In a two-dimensional lattice, one 
can achieve arbitrary displacements of the wave packet 
by suitable modulation of the external field [14]. 

Also excitations (excitons) in idealized linear polymers, 
modeled by beads and springs, with an external field and 
at ultra-cold temperatures have been shown to exhibit 
Bloch oscillations [15]. Similar to the experiments men¬ 
tioned above, one can manipulate the excitons’ motion 
by varying the external field [16]. 

However, all the experiments so far consider lattice¬ 
like underlying potentials. For more complex arrange¬ 
ments one encounters the situation of a (quantum) graph 
[17, 18], where vertices are connected by arms. In order 
to be able to manipulate a wave packet on an arbitrary 
graph, one first has to understand the behavior of wave 
packets at bifurcation points. These points are charac¬ 
teristic for vertices where three (or more) arms meet. We 
call a graph with a single bifurcation a star graph. In the 
following, we will assume that each arm is modeled by a 
one-dimensional periodic potential and that at the vertex 
there is continuity and current conservation, see below. 

As we will show, it is possible to manipulate a wave 
packet by suitable changes of the external fields. Since 
complex graphs can be build iteratively by joining star 
graphs, our results are a first step to precise manipu¬ 
lations of wave packets on these complex graphs. This 
paper is organized as follows. In the next section we will 
recall the problem of a tilted lattice solved in a continu¬ 
ous model. Section III presents the study of the directed 
transport in a one dimensional lattice. In section IV we 


treat such transport in driven quantum star graphs. Fi¬ 
nally, section V presents concluding remarks. 

II. BLOCH OSCILLATIONS IN A 
ONE-DIMENSIONAL LATTICE 

We start by recalling the effect of Bloch oscillations in 
one-dimensional periodic potentials V ( x ). The Hamilto¬ 
nian is given (in units h = 2m = 1) by: 

H =--^-2+ V ( x ) + f x ’ (!) 

where / is the external field strength and V ( x+d ) = V ( x ) 
is the periodic potential with lattice period d. This 
Hamiltonian describes, for instance, an electron in a one¬ 
dimensional crystal in the presence of a constant elec¬ 
tric field [8, 19]. Such a Hamiltonian generates time- 
periodic oscillations, the Bloch oscillations, of an initial 
Gaussian wave packet (GWP). These oscillations have a 
period inversely proportional to the field strength, i.e., 
Tb = 2ivh/df, and a well-defined amplitude L = A//, 
where A is the bandwidth. In the following we choose 
the lattice periodic potential as a cosine potential of the 
form: 

V(x) = Vo cos • ( 2 ) 

Exemplarily, we show in Fig. 1 the Bloch oscillations of 
the solution of the Schrodinger equation with H and an 
initial GWP over two Bloch periods. We chose the pa¬ 
rameters such that our results are comparable to one of 
the discrete model discussed in Ref. [15]. We note that 
the solution of Eq. (1) shows the different bands of the 
Bloch band structure [8] , while in the discrete model only 
the lowest band was considered. However, also for the 
continuous solution, the main fraction of the excitation 
is bounded to the lowest band, see the weak contribution 
of the higher bands with larger velocity in Fig. 1. 
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FIG. 1. Contour plot of the probability density. Bloch oscil¬ 
lations for / = 0.2, d— 1 and Vo = 16.7875. One can see that 
the oscillation is bounded within the interval of the length 
A = A// = 20. 


III. PERIODIC POTENTIALS WITH FIELD 
MODULATIONS 


Now, modulations of the potential or of the external 
field allow to manipulate the motion of the GWP. In 
order to see this, we assume either a time-dependent po¬ 
tential V(x,t) or a time-dependent and spatially homo¬ 
geneous external field F(t). Then, the Schrodinger equa¬ 
tion in the one-dimensional case has the general form 

d 

i ~Q t ^^ x ' t ) = + F (x,t) + F(t)x] V(x,t) . 

( 3 ) 

A simple sinusoidal change in the coupling strength is 
reflected by an oscillating potential amplitude leading to 
tunneling matrix elements of H given by 


V(x, t) = Vo [1 — cl sin ( u)t + <p)\ cos 



( 4 ) 




where a G [0,1]. Fig. 2 shows three cases of the dy¬ 
namics of a GWP centered at xq = 78 with a standard 
deviation a = 6. Each panel displays the motion for a 
constant external field of strength / = 0.2 and an oscil¬ 
lating potential with parameters d = 1, a = 0.85, and 
Vo = 16.7875, but with different phases (p. These param¬ 
eters are such that the results match those of Ref. [15], 
where the oscillating coupling strength was achieved by 
varying distances between the centers of the (discrete) 
potential. 

Experimentally it might be easier to modulate the ex¬ 
ternal field than the potential. Thus, we consider in 
the following only the case of a time-dependent field 
F(t) = f sm(ut + <p ) and a time-independent periodic 
potential V{x) = Vo cos[(27r/d)x]. Note, that the phase 
<p is now also shifted to the field dependence. For peri¬ 
odically modulated external fields, it is known from the 
discrete model that a GWP with an initial width cro can 


FIG. 2. Lattice potential amplitude modulation given by (4): 
Contour plot of probability density for a = 0.85 and uj = df. 
Each plot corresponds to different values of p = 0,7r/4 and 
7t/2. 

broaden with time [12]: 

cr{t) = oo^l + e[Jb(//w)cos[(//w)cos^]/(jg] 2 , (5) 

where Jo(f/w) is the Bessel function of first kind. There¬ 
fore, whenever the Bessel function or the cosine function 
are zero, it is possible to preserve the width of the ini¬ 
tial GWP, leading to a non-dispersing wave packet for a 
GWP with / = 7r/10,c j = 1/5, and <p = 0. In the fol¬ 
lowing we will always use these parameters to preserve 
the width of the GWP (before the bifurcation point). 
Strictly, Eq. (5) is only valid for the lowest band, while 
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contributions from higher bands will in any case result in 
a broadening of the GWP. We confirm the preservation 
of the width of the GWP in the lowest band by numerical 
fitting to a Gaussian with parameterized width. 

Fig. 3 shows the dynamics of the same initial GWP as 
in Fig. 2 with parameters / j= 7t/10, d = 1, uj = 0.2, and 
Vo = 16.7875, for three different phases As for the 
case where the potential is time-dependent, one clearly 
observes the ^-dependence of the dynamics: While for 

= 7r/2 there is no net displacement after one Bloch 
period, one can achieve a significant displacement for f = 
7r/4 and 7r = 0. We note, that for time-dependent fields 
with a phase 0, also the Bloch frequency depends on <fi. 


IV. DIRECTED TRANSPORT IN DRIVEN 
STAR GRAPH 

In this section we solve the problem for a quantum star 
graph. Each of the arms is vested with a tilted lattice 
potential, and they are connected at a common (central) 
vertex. We assign the coordinate x to each arm, which 
indicates the position along the arm; x takes the value 0 
at the common vertex. 

Thus, we consider the time-dependent Schrodinger 
equation for a star graph with N arms (h m 2 m = 1): 

(6) 

Vt, 0 < x < Lj , j = 1,..., N 

where 

Vo cos . 

As mentioned in the previous section, the time- 
dependence of the field is given by Fj(t) = fj sm(ut+<fij), 
being in this case also arm-dependent, and the parame¬ 
ters fj and (fj are used to manipulate the GWP transition 
through the central vertex. The imposed conditions 

' *i(L u t) = 4/ 2 (L 2 , t) = ... = 9 N (L N ,t) = 0 , 

= * 2 ( 0 , t) = ... = 4 / tv ( 0 , t) , , . 

\ N (7) 

E £*j(x,t )\ x=0 = 0 . 

I 3 = 1 

imply Dirichlet boundary conditions at the non- 
connected edges and continuity and current conservation 
at the vertex. 

The solution of Eqs. (6)-(7) can be written in terms 
of the complete set of eigenfunctions 'ipj^ n (x) of the 
potential-free star graph: 

= ^2c n (t)'tl) j: n(x), j = l,...,N. (8) 
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FIG. 3. The external field modulation given by F(t) — 
f sm(ojt+(/>): Contour plot of probability density for / = 7r/10 
and uj ass 0.2. Each plot corresponds to different values of 
<f = 0,7r/4 and 7r/2. 


For the star graph with N arms the eigenfunctions 
ifj, n (x) of the stationary Schrodinger equation (in units 
k = 2m = 1): 

d 2 

0 < x < Lj, j = 1, (9) 

with the same boundary conditions as in Eq. (7) have the 
form 


B„ 


sin (k n Lj) 


sin [k n (Lj -x)], 


( 10 ) 


where time-dependent coefficients C n (t) are to be found. 


ipj,n(x) 
































































































































4 


where 


account the orthonormality of the eigenfunctions, 


B 


n 


- - 1/2 

T J [Lj + sin (2 k n Lj)\ sin -2 (k n Lj)/2 (11) 

j 



j j ,ndx — $mm 


(13) 


we get the system of ODE with respect to the coefficients: 


are the normalization coefficients, and the k n fulfill the 
equation: 


N 

cot(k n Lj) = 0. (12) 

3 = i 


Eq. (12) can be solved numerically. In order to avoid non¬ 
generic degeneracies, the lengths of arms Lj are chosen 
to be rationally independent [17]. 

Inserting the expansion (8) into Eq. (6) and taking into 


— k n C n + ^ ^ M nrn C m . (14) 

m 

The matrix M can be written in terms of the contri¬ 
butions arising from the lattice potential (/^) and the 
driving external field (/( F )): 

M = / (y ) +/ (F) , 

the elements of which are given in Appendix A. 

As the system of ODE (14) has an infinite number 
of equations, solving it numerically requires to limit this 
number. Then, to preserve the required accuracy one has 
to control the norm conservation condition 

J2\Cn(t)\ 2 = l, Vi. (15) 



FIG. 4. Contour plot of the probability density for the three arm star graph, where driving external field is given by Fj(t) = 
fj sin (ujt + 4>j ) with fi = —f <2 — —fz = 7r/10, cj = 0.2. The three columns correspond to the arms of the star graph. The 
^-coordinate of the first arm is reversed. 


We now turn to the results for a GWP on a star graph. 
For clarity, we restrict ourselves to a star graph with three 
arms. Fig. 4 shows the dynamics of a GWP initially 
located on the first arm, at xo = 22, see the leftmost 
panel in Fig. 4. In this first example, the external field 


strengths Fj(t) = fj sin (ut + <fj) are of the same magni¬ 
tude / for all the arms and chosen to be /i = = — fs, 

such that the field points globally in one direction. As 
for the case of Bloch oscillations on a single arm, the 
GWP moves towards the vertex (at position 0). At the 
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FIG. 5. The time-dependence of the partial norms corre¬ 
sponding to Fig. 4. 


vertex the GWP bifurcates and is partly transmitted to 


the other two arms and partly reflected back into the 
initial arm. As the driving external field for the second 
and third arms has the same parameters, the dynamics 
of the split GWP is identical for these two arms. This 
can also be seen in Fig. 5, where the time-dependence of 
the partial norms 

Pj(t) = J\^j(x,t)\ 2 dx, j = 1,2,..,7V 
0 

are shown. We have confirmed that the total norm is 
= 1- On the basis of previous results for 
quantum graphs with the Neumann boundary condition 
[17], we can calculate the fraction of the GWP which is 
reflected back into the initial arm as being 1/9. The other 
fraction of 8/9 of the initial GWP is split into the two 
identical parts. Due to the reflection and transmission 
at the vertex the shape of the GWP in the two arms 
is (slightly) distorted, but with a width which is still 
roughly preserved over time. 



FIG. 6. Contour plot of the probability density for the three arm star graph, where driving external field is given by Fj{t) — 
fj sin (wt + 4>j ) with fi — —/ 2 = —f$ — 7r/10,u; = 0.2. The three columns correspond to the arms of the star graph. The 
^-coordinate of the first arm is reversed. 


The situation changes when we change the phase of 
the second arm to </>2 = 7t/ 2, see Fig. 6. While the trans¬ 
mission into the third arm persists, although the wave 
packet becomes more distorted, the transmission in the 
second arm is nearly prohibited. This effect can be quan¬ 
tified by calculating the time-dependence of the partial 
norms, which are shown in Fig. 7. After a transient time 
until most reflection and transmission effects have taken 
place, the partial norm in arm three saturated around a 
value of 70%, while the partial norms of the other two 


arms saturate around values between 10 and 20 percent. 
Even though the wave packet in the third arm gets dis¬ 
torted it still retains a rather constant width. Therefore, 
by manipulating the phase of the external field in one 
of the arms, one is able to “block” the transmission into 
this arm and consequently can manipulate the direction 
of the initial GWP. 

For our setup of initial parameters, the choice (j )2 = 7r/2 
turns out to be the most efficient one to realize the idea 
of “blocking” the wave packet propagation to the second 
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FIG. 7. The time-dependence of the partial norms corre¬ 
sponding to Fig. 6. 



FIG. 8. The phase-dependence of the partial norms at time 
t ss 94.25. The external field and initial parameters corre¬ 
spond to Figs. 4 and 6. 


arm. To show this, we calculated the partial norms for 
the three arms depending on the phase 02- Fig. 8 shows 
the value of the partial norms at times where the norms 
have saturated, i.e., we have chose that time at which 
we stop our numerical calculations, t = 94.25. At the 
value of (fi 2 = 7t/2, we observe the largest difference in 


the partial norms. 

V. CONCLUSION 

In this paper, we have investigated the possibilities of 
directing a Gaussian wave packet on a star graph with 
three arms. This can be viewed as a paradigmatic build¬ 
ing block for complex networks, which can be composed 
by iteratively joining star graphs. Assuming that each 
arm of the graph is equipped with a periodic potential 
and additionally (independent) possibly time-dependent 
external fields act on each arm, our numerical results 
show that it is possible to direct a wave packet with high 
probability from one arm into one of the other arms while 
blocking the transmission into the third arm. This ef¬ 
fect is best achieved by assuming time-periodic external 
fields with independent phase shifts. The phase shifts 
allow for a parametric tuning of the “blocking” effect. 
Our results further indicate that the shape of the wave 
packets remains almost intact. We believe that such a 
scenario can be realized by state-of-the-art experiments 
with, say, ultra-cold atoms in optical lattices. One can 
also imagine that after transmission of the wave packet 
into one arm, the shape of the wave packet can be re¬ 
established, see, e.g., [20]. The fact that there is not 
perfect transfer into the desired arm can be explained by 
fundamental reflection and transmission rules at the bi¬ 
furcation point, see also [17]. This study being a proof 
of concept, we will extend our analysis to more complex 
situations where several star graphs are joined together 
forming more complex networks, such as, e.g., T-fractals 
[21] or dendrimers [22]. 
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Appendix A: Calculation of matrices and I^ 


The elements of the matrices 1^ and /^ have integral forms, which can be solved analytically. Denoting the 
frequency of the periodical lattice potential by ujd = 2it/d one gets the following: 
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N L J 


Lm = ^0 Ylj Iplm C0S {U d x)lpj, n dx (Al) 

"sin (u d Lj) + sin (k n - k m )Lj sin ( w d Lj ) - sin (k n - k m )Lj 


i- 1 o 
N 


= V 0 J2 


B n B r , 


4 sin ( h n Lj ) sin (^k^Lj^j _ H - k n /c m uj^ k n H - 

sin (ujdLj) + sin (fc n + k m )Lj sin (oo^Lj) — sin (fc n + k m )Lj 


+ k n + 


kn k 75 




AT J AT 

E F iW / ^tm x ^i,ndx = J2 
j =l n 


Fj^BnBn 


—[ sin (fenLj) sin ( k m Lj) 


where the diagonal elements of matrix ^4 (i.e. for m = n) are: 

1 — cos (2k n Lj) 

Ann = Y 4^1 

and the off-diagonal elements (m ^ n) are: 

1 — COS 2 (fc n — k^Lj 1 — COS 2 (fc n + k m )Lj 


A — 


(&n k m )2 


(^n 4“ k m )2 


(A2) 
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